Ferromagnetic spin fluctuation induced superconductivity in Sr 2 Ru04 
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We propose a quantitative model for triplet superconductivity in Sr2Ru04 based on first principles 
calculations for the electronic structure and magnetic susceptibility. The superconductivity is due 
to ferromagnetic spin fluctuations, that are strong at small wave vectors. The calculated effective 
mass renormalization, renormalized susceptibility, and superconducting critical temperature are all 
in good agreement with experiment. The order parameters is of comparable magnitude on all three 
sheets of the Fermi surface. 
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The layered ruthenate, Sr2Ru04 has attracted con- 
siderable recent interest. It is structurally similar to 
the first cuprate superconductor, (La,Sr)2Cu04, is near 
a magnetic instability (Sr^Cai-aRuOg and Sr 2 RuYC>6 
are ferro- and antiferromagnetic, respectively), and was 
thought to be strongly correlated. However, closer ex- 
amination reveals more and more differences from the 
cuprates. It was noted that SrRuC>3 is ferromagnetic 
(FM) so it was conjectured that Sr 2 Ru04 must be close 
to a FM instability as well KJ. This has recently been 
corroborated by detailed microscopic calculations of mag- 
netic properties of ruthenates . Ferromagnetic fluctua- 
tions disfavor both s— and d— wave superconductivity, so 
it was suggested |||| that superconductivity in Sr 2 RuC>4 
must be triplet (p-wave) , thus very different from the high 
T c cuprates. The idea of the strongly correlated electrons 
in Sr 2 Ru04 is mostly based on the apparent disagreement 
of Angular Resolved Photoemission (ARPES) measure- 
ments of the Fermi surface Q] with results of LDA cal- 
culations HD- This argument is, however, questionable, 
because ARPES measurements disagree with de Haas- 
van Alphen (dHvA) experiments |tJ. The notorious fail- 
ure of the LSDA to describe properly the magnetism of 
undoped cuprates does not occur in ruthenates . Thus, 
the case for the strong correlations in ruthenates is ques- 
tionable. 

The following challenges should be met by a quantita- 
tive theory of the electronic states and superconductivity 
in Sr 2 Ru04: (i) reconciliation of the (well reproducible) 
ARPES results with the dHvA measurements. This may 
also relate to whether or not there are strong correla- 
tion effects, (ii) the mechanism for superconductivity 
and how it is related with the large mass renormaliza- 
tion (of a factor of 3-4). In this Letter we address both 
of these issues. 

The valence bands of Sr 2 RuC>4 are formed by the three 
tig Ru orbitals, xy, yz, and zx. These are hybridized with 
the in-plane oxygen and, to a lesser extent, with the api- 
cal oxygen p-states. The bare oxygen p levels are 
well (~ 2 eV) removed from Ep, so the effect of the O 
p orbital is chiefly renormalization of the ionic i 2g lev- 
els, and assisting in the d — d hopping. The LDA band 



structure can be reasonably well described in the vicin- 
ity of the Fermi level as three mutually non-hybridizing 



tight-binding bands: e xy (k) = 
cos ak y ) + 4t' ddv cos ak x cos ak y , 



E + 2t dd7r (cosak x 
and e^*y x (k) = E 



2t ddAZaf} + 8 ^ cos cos cos 1h with the Pa- 



rameter (E -E F ,t ddv ,t' dd7r ,t ± ) being (-0.4, 0.4, -0.12,0) 
and (-0.3, 0.25, 0, -0.025) eV for the xy and xz, yz 
bands, respectively, for the bands of Ref. [^J. With near- 
est neighbors only, this model yields one nearly circular 
cylindrical electronic sheet (7) of the Fermi surface (FS) 
and four crossing planes (quasi-lD FS). The weak xz — yz 
hybridization reconnects these planes to form two tetrag- 
onal prisms, a hole one [a) and an electron one (/?), as in 
Fig.0. ARPES gives a FS of different topology: the van 
Hove singularity at k= (ir/a,0), which appears slightly 
above the Fermi level in the calculations (« 60 meV in 
the LDA calculations and ss 50 meV with the gradient 
correction [g included), is seen below it in photoemission 
experiments. This reconnects the surface 7 and makes 
it hole-like instead of electron-like. The total electron 
count in the ARPES FS is still 4, indicating stoichiom- 
etry of the samples. Important consequences were as- 
cribed to the fact that the van Hove singularity is situated 
in the same place as in the cuprates. The main differ- 
ence between the LDA and the ARPES Fermi surfaces is 
that the latter corresponds to a larger _E ({^}z) — E(xy) 
(the levels are higher because of an additional hy- 

bridization with the apical (02) oxygen P^y orbitals). 
Importantly, the calculations imply strong Stoner renor- 
malization. The Stoner factor /, calculated as described 
in Ref. @], is 0.43 eV, and AT(0) = 2.06 eV" 1 §. This 
yields a Stoner renormalization 1/(1 — IN) — 9, some- 
what larger than deduced from experimental susceptibil- 
ity x/Xband = 7.3. To fit the experiment, / should be 
Icxp = 0.42 eV. Note that the experiment leaves no room 
for any renormalization of \ beyond the Stoner one. 

However, the topology of the ARPES FS disagrees with 
that from dHvA experiments. The latter yields three 
cross-sections, which sum up to 4 electrons/cell with ex- 
cellent accuracy only if the surface 7 is electron-like. The 
LDA calculated a, (3, and 7 areas deviate from the dHvA 
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experiment by only -2%, -3% and 5% of the Brillouin zone 
area, respectively, and an exact match can be achieved 
by very slight shifts of the bands a, (3, and 7 by 5, -4, and 
-3 mRy, respectively. Such agreement is generally con- 
sidered very good even in simple metals, and the small 
mismatch (which does not change the FS topology) is 
likely due to some underestimation in LDA calculations 
of the tiny xz — yz hybridization. Unlike dHvA, which 
probes the bulk, ARPES probes essentially first surface 
Ru-0 layer. The cleavage plane in Si~2Ru04 is likely as- 
sociated with the rocksalt layers, leaving the Ru-02 bond 
dangling or otherwise strongly perturbed. As such, this 
bond is likely to be contracted compared to the bulk, and 
the electronic structure of the surface RuC>2 layer differs 
from bulk. The main effect of such a surface relaxation 
is expected to be a strong modification of the Ru({^}z)- 

02 ({y}) hopping. In a linear approximation, this can be 
estimated from bulk calculations with the Ru-02 bond 
length reduced by a half of the supposed surface contrac- 
tion of this bond. We performed such calculations for 
Sr2Ru04 with the 02 shifted by 0.1 A and found that 
the energy distance between the Fermi level and the van 
Hove singularity was reduced by 30 meV. Thus, the sur- 
face relaxation of the Ru-02 bond that would bring our 
LAPW calculations in agreement with the ARPES mea- 
surements would be less than 0.4 A (probably, closer to 
0.3 A, due to non-linearity Although the actual sur- 

face relaxation for Sr2Ru04 is not known, the change in 
the observed electronic structure due to the Ru-02 bond 
relaxation is in the right direction and of the right order 
of magnitude compared to the observed ARPES FS. We 
conclude that the LDA and dHvA yield the bulk elec- 
tronic structure of Sr2Ru04. The differences in ARPES 
presumably reflect the surface structure. 

This said, we recall that the mass renormalization 
(1+A) found in dHvA experiments and from the spe- 
cific heat is unusually large: for the xy (7) sheet 
it is 4, and 3.3 for the two other sheets. Materials 
with an electron-phonon coupling constant of the order 
of 2.5 are known, but if it were so large in Sr 2 Ru04, 
with its high phonon frequencies, the superconducting 
T c would be much higher than 1.5 K. This paradox is 
naturally resolved in the framework of the conjecture 
jl],H that Sr 2 Ru04 has strong electron-paramagnon cou- 
pling, and may even be a p-wave superconductor, which 
is also in accord with recent experiments showing anoma- 
lously strong dependence of T c on residual resistivity jl2) . 
In such a case, two different coupling constants appear: 
A™ which controls the mass renormalization is the aver- 
age of the electron-paramagnon interaction over the FS, 
while A™ 1 which determines the p-wave transition temper- 
ature is the I = 1 angular component of this interaction. 
Importantly, this holds for any boson-mediated interac- 
tion, including electron-phonon, so that A = A™ + Aq, 
Ai = A" 1 + A^. For s-pairing although the mass renor- 



malization is controlled by A™ + Aq, the superconducting 
coupling constant is Aq — A™. 

The situation with Sr2Ru04 is further complicated by 
the fact that there are three different sheets of the FS 
and the order parameters on all sheets should be deter- 
mined simultaneously. Without a quantitative numerical 
estimate it is impossible to assess whether or not the 
triplet pairing hypothesis of Refs. [QJ|] can be reconciled 
with the body of experimental facts. Fortunately, the 
LSDA calculations provide the necessary information for 
a quantitative analysis. 

The most important (and most uncertain) part of such 
an analysis is the interaction responsible for pairing and 
for the mass renormalization. This was not specified in 
the previous works, but we conjecture that it is the ex- 
change of paramagnons. Such an interaction in metals 
was studied with respect to possible superconductivity 
in Pd in the late 70-ties (see, e.g., Jl3|,[l4[] ) , and later in 
connection with heavy fermions. The parallel-spin inter- 
action, relevant for triplet pairing is given in the RPA by 
the sum of the diagrams with odd numbers of loops, 



V(q = k - k') 



I 2 (q)Xo(q) 



(1) 



The mass renormalization is not as easy to define. 
Besides the parallel-spin interaction (|lj), there is the 
antiparallel-spin interaction, given in the same approx- 
imation by the sum of the chain diagrams with even 
numbers of loops, plus ladder diagrams |l^,|l^]. In case 
of contact interaction, the total interaction is three times 
stronger than the interaction in the parallel-spin chan- 
nel only. It was pointed out though, that there 
is no good physical reason to single out any particular 
class of diagrams. It was found that including all three 
classes above leads to systematic overestimation of mass 
renormalization by a factor of 2 to 3 (R]|l(| • Our case is 
further complicated because unlike the electron-phonon 
interaction, the electron-electron (and, correspondingly, 
the electron-paramagnon) interaction is already included 
in some average way in the LSDA band structure. Thus, 
the electron-paramagnon mass renormalization is to some 
extent included in the LDA mass as well. 

Despite all this difficulties, one can get an idea about 
the size of the electron-paramagnon mass renormaliza- 
tion by making calculations with the parallel-spin inter- 
action (Q) only; one may consider that as a lower bound 
for the total spin-fluctuation induced renormalization. 
The mass renormalization then is computed in the same 
way as the electron-phonon renormalization, i.e, by tak- 
ing the average of V(q) of Eq.(Q) over the FS. One has 
to remember, though, that there are other effects beyond 
LDA, apart from the one that we calculate, which may 
further increase the observable mass. 

The triplet pairing constant is calculated by averag- 
ing V(q) with the functions reflecting the k-dependence 
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of the direction of the (vector) order parameter, in the 
simplest case with k • k'/fc/c'. A common approximation, 
which we use here (although it may be not as good in 
Sr 2 Ru0 4 as in Pd) is to take xo(l) = Xo(0) = N. The q 
dependence of / cannot, however, be neglected and has 
to be specified. Essentially, it tells us how much the FM 
state is favored over AFM states p7| . As discussed in 
Ref. 0, what favors ferromagnetism over antiferromag- 
netism in ruthenates is the oxygen contribution to the 
Stoner factor. This is determined from the band struc- 
ture calculations as follows: Atomic Stoner factors for 
Ru and O ions are calculated in a standard way and 
are Ir u « 0.7 eV, Io ~ 1-6 eV. The total Stoner fac- 
tor for the compound is / = Ir u v% u + 21oVq, where 
vr u and vo are partial densities of states at Ep of Ru 
and in-plane oxygen; the contribution of the apical oxy- 
gen is negligible. For AFM ordering, the second term 
in the expression for / falls out, because oxygen is non- 
magnetic by symmetry. We found the AFM Stoner fac- 
tor / for Sr 2 RuC>4 to be smaller than FM one by 14% 
(oxygen contribution AI = 0.06 eV). A g-dependence 
that reflects this effect is I(q) — 1/(1 + b 2 q 2 ), where 
b 2 = 0.5(a/7r) 2 A//(/ - AI) w 0.08(a/7r) 2 . In the fol- 
lowing we use this I(q) together with EqjlJ 

Let us now make link to the real FS. In Refs. [p]Jl8| 
the maximum full cylindrical symmetry was assumed 
for all three FS's. This approximation completely ne- 
glects the quasi-lD character of the xz and yz bands 
and cannot be used for quantitative purposes. Instead, 
we retain the cylindrical approximation for the xy FS 
7 and use the ID approximation for the xz and yz 
FS's. Then we have three 2D Fermi lines: 7, a cir- 
cle with the radius g ~ 0.97r/a, £, two lines parallel to 
x at w ±2<7 / 3 from the T point, and £, the two cor- 
responding lines parallel t o y (FigJ^). Using the stan- 
dard multiband technique |Q we now introduce the cou- 
pling matrix A£ = N~ l Y^uk' S(eia)S(evj)V(k. - k') = 
NviVj < V(k — k') >ij, where (i,j) can be 7, £, or £, 
and Vi = Ni/N (from our band structure vf~ ( rj w F{,,o 
and v 1 = 0.44, = 0.28). Then the mass renormaliza- 
tion in band i is defined as Af = ^ . Ay. The average 
mass renormalization is X s — V Af„- . 

Using this model, we arrive at A*^ = 0.35, A|t = 0.32, 



A 



^ = 0.16, A| c = 0.03. This gives A* 



(Af 



1. 



2A^)/^ = 1.5, A| = (A| ? + A^ 
A s = 1.7, to be compared with experimental dHvA val- 
ues of 3, 2.3, and 3, respectively. The difference may 
be due to an electron-phonon coupling of the order of 
1 and/or antiparallcl spin fluctuations, neglected in our 
calculations. 

Let us now return to the question of the p-wave super- 
conductivity. The theory for a cylindrical FS is presented 
exceedingly well by Sigrist et al (l],^(],|l8| and need not be 
repeated here. The only difference for a FS of arbitrary 
shape is that instead of the k- vector components, we have 



to use Allen's FS harmonics |19| . So, there are four pos- 
sible unitary planar states, all degenerate if spin-orbit is 
neglected. Let us consider, for instance, the A\ u state: 



d k = d — , 



(2) 



where Vk is the Fermi velocity. This state has, generally 
speaking, a finite superconducting gap, and thus zero 
density of states at the Fermi energy below T c , in contrast 
with the experiment |ll[ . The same holds for three other 
states, degenerate with the one of Eq.(||). Nonunitary 
linear combinations of the states which are gapless are 
also possible. These have, however, generally speaking, 
lower pairing energy and should not occur. 

We now calculate the transition temperature within 
our spin-fluctuation model. Similar to Agterberg et al 
fTsfl , we consider the superconducting state with the or- 
der parameter d which is constant for each of the three 
FS sheets, but differing between the sheets. We have 
to calculate the matrix A p - = Nv. k Vj < V{k — k')(dj c ■ 

dk')/(^k^k') where i and j label the three bands, 
and find the maximum eigenvalue of the matrix f,- A?. 
|]l9| . The corresponding eigenvector defines the coeffi- 
cient a and thus the relative magnitude of the order pa- 
rameter in bands 7 and (£, £). By symmetry, the pairing 
matrix looks like 

AL A^ 
AL 
A p A p 

(If we had used instead of £,£ nomenclature the a, f3 
one, as in Ref. |T§|], this symmetry would not hold). 
Numerical calculations give A^ 7 = 0.16, A?* = 0.075, 
and A^£ = 0.025. The maximum eigenvalue of the corre- 
sponding coupling matrix is 0.43, and the corresponding 
superconducting state is 0.857 + 0.38£ + 0.38C- 

Let us now estimate the transition temperature. Us- 
ing the characteristic paramagnon energy u s f ~ (A -1 — 
I)/4 « 160 K from our calculations, as the cut-off fre- 
quency, and the Allen-Dynes formula for strong coupling 
(although A = 0.43 is relatively weak, the relevant num- 
ber is the renormalization parameter X s ks 1.7), we obtain 



T, 



(N- 1 - J)/4 
1* 



exp[-(l + A s )/AL f ] = 0.25K. (3) 



Again, as in the case of mass renormalization, there is 
some room for the electron-phonon coupling as well. 

One of the key problems, as discussed in Refs. p0| , |l8| |, 
is the residual electronic specific heat M], which re- 
mains at about 50% of its normal value in the super- 
conducting regime. This led Agterberg et al [M to pos- 
tulate a pairing matrix which yields a vanishing gap for 
the 7 band. This, however, does not square with the 
quantitative estimate presented in this Letter. An ear- 
lier assumption |^|| was that the excess pairing en- 
ergy that forbids nonunitary combination of the order 
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T e dlog N 
dE F 



parameters (g) may be overcome by additional magnetic 
(Stoner) energy in a nonunitary state. The requirements 
are strong Stoner renormalization ( supported by the 
calculations) and strong particle hole asymmetry pD. 
However, a quantitative estimate according to Ref. |2l| ] 
shows that the effect is by far too weak. The criterion 

jztn log ~ 10 -5 , while it should be of 
the order 1 for the nonunitary state to stabilize. So, the 
problem of the residual electronic specific heat remains 
open, although it cannot be excluded that it results from 
sample inhomogeneity and is extrinsic to superconduc- 
tivity. 

To summarize, we have presented first principles cal- 
culations which indicate that (1) conventional LDA cal- 
culations give a correct description of the bulk electronic 
structure of SrgRuO^ as well as of its renormalization 
due to the Stoner exchange interaction; (2) the difference 
between the bulk and the surface electronic structure (as 
measured by ARPES) can be explained by the surface ef- 
fect; (3) interactions due to exchange of FM spin fluctua- 
tions, as calculated from the LDA band structure, are suf- 
ficiently strong to explain both the mass renormalization 
and superconducting critical temperature. We would like 
to emphasize the main approximations used: (a) Neglect 
of the q-dependence of xo > (b) neglect of strong coupling 
effects beyond the Allen-Dynes formula (note that strong 
coupling effect will lead to a final density of states N(0) 
at all T and thus to nonexponential T-dependencies 
of specific heat and like quantities), and (c) neglect of cor- 
relation effects in the mass renormalization beyond the 
parallel-spin paramagnon induced interaction. In princi- 
ple, it is clear how to improve the first two items, while 
the last issue lacks full theory and cannot be easily dealt 
with. 

We acknowledge helpful discussions with W. Pickett, 
R. Rudd, and D. Hess. Work at NRL is supported by 
the ONR. Computations were performed at the DoD 
HPCMO NAVO Center. 



[1] T.M. Rice and M. Sigrist, J. Phys. Condens. Matter 7, 
L643 (1995). 

[2] I.I. Mazin and D. J. Singh, (Phys. Rev. B, in press). 
[3] K. Machida, M. Ozaki, and T. Ohmi, J. Phys. Soc. Jpn. 

65, 3720 (1996). 
[4] T. Yokoya et al, Phys. Rev. Lett, 76, 3009 (1996); D.H. 

Lu et al, ibid, 76, 4845 (1996). 
[5] D.J. Singh, Phys. Rev. B 52, 1358 (1995). 
[6] T. Oguchi, Phys. Rev.B 51, 1385 (1995). 
[7] A.P. Mackenzie et al, Phys. Rev. Lett. 76, 3786 (1996); 

ibid, 78, 2271 (1997). 
[8] J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 

(1992). 



[9] Y. Maeno et al, Nature 372, 532 (1994). 

[10] Since oxygen p bands are well below Ef, the effect is 
proportional to t 2 , where t is the pdn Ru-02 hopping, 
for which the canonical scaling is t oc d~ 4 . 

[11] Y. Maeno et al, J. Low Temp. Phys. 105, 1577 (1997). 

[12] Y. Maeno, to be published. 

[13] D. Fay and J. Appel, Phys. Rev. B 22, 3173 (1980), and 

references therein. 
[14] P.B. Allen and B. Mitrovic, Solid State Phys., 37, 1 

(1982). 

[15] S. Doniach and S. Engelsberg, Phys. Rev. Lett. 17, 750 
(1966). 

[16] K. Levin and O.T. Vails, Phys. Rev. B 17, 191 (1978). 

[17] We do not discuss here the interplanar magnetic interac- 
tions, which, contrary to the intraplanar one, are weak 
and probably antiferromagnetic Q. 

[18] DF. Agterberg, T.M. Rice, and M. Sigrist, Phys. Rev. 
Lett. 78, 3374 (1997). 

[19] W.H. Butler and P.B. Allen, in: Superconductivity in d- 
and f -metals, ed. by D.H.Douglass (Plenum, N.Y.,1976). 

[20] M. Sigrist and M.E. Zhitomirsky, J. Phys. Soc. Jpn. 65, 
3452 (1996). 

[21] T. Sugiyama and T. Ohmi, J. Phys. Soc. Jpn. 64, 2746 
(1995). 



X , . . , X 




X ' " ~ ' X 

FIG. 1. LAPW Fermi surface. The origin is at P The 
thickness of the lines is inversely proportional to the Fermi ve- 
locity (the inner and the outer contours are Ef±2 mRy). The 
model Fermi surfaces £ and C, are shown as straight lines. The 
model Fermi surface 7 is within the 2 mRy window around 
the actual surface and thus not shown. 
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